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Abstract 

An  ultracapacitor  model  with  automatic  order  selection  of  complexity  and  automatic  scaling  of  capacity  is  created  in  Virtual  Test 
Bed  platform  for  complex  system  simulation  and  prototyping.  The  order  selection  is  based  upon  the  variable  simulation  time  step  for 
appropriate  level  of  details  of  the  model,  and  automatically  executed  by  the  model  itself,  therefore  it  applies  to  simulation  for  both  high  and 
low  frequencies,  or  for  both  fast  and  slow  transients,  or  even  for  multi-time  scale  resolutions.  A  simple  trapezoidal  algorithm,  rather  than 
the  general  stiff  algorithm,  which  is  complicated  in  implementation,  is  used  for  numerical  integration.  For  a  given  time  step,  the  truncation 
error  is  controlled  to  be  negligible,  so  that  numerical  stability  and  accuracy  are  ensured  for  the  given  order  of  model  complexity.  The  model 
is  validated  by  comparing  the  simulation  results  to  the  experimental  ones.  The  application  example  of  the  ultracapacitor  model  in  complex 
system  simulation  is  also  presented. 
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1.  Introduction 

An  ultracapacitor  or  an  electrochemical  capacitor,  or  a 
double  layer  capacitor  [1]  can  have  much  greater  charge 
or  energy  storage  than  a  conventional  capacitor  because  of 
its  highly-amplified  surface  area  of  micropore-structured 
electrodes.  For  the  same  reason,  an  ultracapacitor  exhibits 
operating  characteristics  that  are  distinct  from  those  of  a 
conventional  capacitor.  In  particular,  its  volt-second  relation 
is  nonlinear,  and  its  impedance  is  a  complicated  func¬ 
tion  of  frequencies.  The  peculiarities  of  the  ultracapacitor 
characteristics  arise  from  the  charge  (anions  and  cations) 
transport  process  in  porous  electrodes  (activated  carbon 
or  other  carbon-based  pore  materials),  which  is  subject  to 
limitations  from  both  mass  transfer  and  Faradaic  process 
that  earned  it  the  name  of  the  electrochemical  capacitor. 
Depending  upon  the  constructions  of  porous  electrodes,  the 
characteristics  of  ultracapacitors  can  vary  significantly  due 
to  different  charge  and  current  distributions  determined  by 
the  pore  structures  and  materials  of  the  electrodes. 

Many  models  for  ultracapacitors  were  developed,  among 
which  is  a  widely  used  one  that  was  based  on  the  porous 
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electrode  theory  [2-4].  The  model  solves  diffusion  equa¬ 
tions  with  respect  to  the  charge  species  transport  in  the 
matrix  phase  of  electronically  conductive  pores  and  an 
ionically  conductive  electrolyte.  The  time-domain  solution 
of  the  model  under  a  galvanostatic  discharge  is  a  combi¬ 
nation  of  a  series  of  transients  characterized  by  different 
time  constants  determined  by  parameters  of  micropore 
structure  and  material  properties.  From  the  stand  point  of 
modeling  and  simulation,  the  significance  of  the  model 
based  on  the  porous  electrode  theory  is  that  it  reveals  how 
the  operating  characteristics  of  ultracapacitors  are  affected 
by  charge  transport  process  for  given  electrode  materials 
and  structures.  In  particular,  an  equivalent  circuit  based  on 
the  transients  can  be  constructed  to  synthesis  the  model 
characteristics,  with  each  transient  being  specified  by  a 
time  constant  generated  from  a  resistor-capacitor  branch. 
Specifically,  a  cascaded  multi-stage  resistor-capacitor  lad¬ 
der  model  can  be  derived  for  the  ultracapacitor  model. 
Depending  upon  the  accuracy  requirement,  the  transient 
terms  can  be  truncated  to  an  appropriate  number  so  that  the 
computation  effort  can  be  minimized  for  a  given  accuracy. 
Correspondingly,  the  number  of  stages  in  the  ladder  model 
can  be  truncated  for  the  same  reason. 

The  ladder  model  derived  by  Miller  et  al.  [5,6]  was  based 
on  the  impedance  spectroscopic  measurement.  The  values 
for  the  resistors  and  capacitors  and  the  ladder  network  can  be 
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viewed  as  a  circuit-parameterized  charge  transport  in  porous 
electrodes  of  ultracapacitors.  In  particular,  the  five- stage  lad¬ 
der  model  has  a  sufficient  accuracy  for  most  applications,  ac¬ 
counting  for  the  frequency  range  up  to  10  kHz,  or  a  transient 
of  a  time  constant  about  100  jjls.  However,  the  five- stage  lad¬ 
der  model  is  often  inconvenient  in  simulation  if  the  system 
characterization  concerns  low  frequencies  or  a  large  time 
constant  only,  since  it  simply  takes  too  much  time  to  calcu¬ 
late  unnecessary  details.  In  this  case,  a  lower  order  or  fewer 
stages  in  the  ladder  are  needed.  Of  course,  a  separate  model 
can  always  be  made  to  serve  the  purpose,  but  our  model 
presented  in  the  paper  is  aimed  to  solve  the  problems  con¬ 
cerning  multi-time  scale  resolutions  with  a  single  model. 

The  model  presented  here  is  based  on  the  five-stage  lad¬ 
der  model  for  the  Maxwell  PC  100  module,  and  it  is  built  in 
the  Virtual  Test  Bed  [7]  computational  environment  for  com¬ 
plex  engineering  system  simulation  and  prototyping.  The 
primary  features  of  the  model  include:  (1)  automatic  selec¬ 
tion  of  the  order  or  number  of  stages  of  the  ladder  network 
according  to  the  simulation  time  step;  (2)  scaled  values  for 
resistors  and  capacitors  in  the  network  when  the  network  or¬ 
der  advances  or  reduces  in  order  to  minimize  the  error;  (3) 
automatic  selection  of  the  network  order  according  to  the 
frequency  input  by  users  and  recommendation  of  appropri¬ 
ate  simulation  time  step;  and  (4)  total  capacity  scaling  based 
on  the  100  F  module. 

In  the  following  section,  we  will  describe  the  model  first 
with  respect  to  the  numerical  accuracy  and  stability,  the  re¬ 
lations  of  the  order  to  the  time  step  size,  and  the  auto-order 
selection  method.  Then  in  Section  3,  the  simulation  of  the 
ultracapacitor  in  an  electric  hybrid  vehicle  is  illustrated.  The 
conclusions  are  given  in  the  final  section. 

2.  Model  description 

The  equivalent  circuit  of  the  multi-stage  ladder  model 
[5,6,8]  for  ultracapacitor  is  shown  in  Fig.  1.  The  reasons  we 
use  this  circuit  are:  (1)  it  physically  mimics  the  distributed 
nature  of  the  ultracapacitor,  and  its  performance  fits  the  ex¬ 
perimental  data  very  well  within  a  wide  range  of  frequencies; 
(2)  it  can  be  easily  combined  with  various  loads  and  used 
to  find  analytical  or  numerical  solution;  and  (3)  the  network 
is  easy  to  adopt  order  reduction  method  as  the  energy  stor¬ 
age  components  in  the  multi-stage  ladder  from  left  to  right 
in  turn  represent  the  complex  transient  processes  from  fast 
speed  to  slow  speed  in  ultracapacitors.  Notice  that  the  model 


Fig.  2.  Comparison  of  the  model  predictions  of  impedance  with  the 
measurements.  The  simulation  results  are  indicated  by  lines  and  the 
measurements  are  by  squares  and  circles. 

contains  actually  six  energy- storage  elements  (including  the 
stray  inductance),  therefore  yielding  six  state  variables.  The 
corresponding  differential  equation  will  be  of  six  orders. 
However,  for  the  convenience  of  discussing  order  reduction, 
we  will  use  the  word  order  and  the  stage  interchangeably. 
For  example,  order  reduction  means  reducing  the  number  of 
stages.  In  Fig.  2,  the  comparison  of  the  impedance  obtained 
from  the  equivalent  ladder  model  (as  in  Fig.  1)  and  from  the 
measurements  [8]  is  shown  for  the  Maxwell  100F  ultraca¬ 
pacitor.  The  parameters  for  the  model,  as  given  by  Table  1 , 
were  obtained  by  using  nonlinear  least- squares  (CNLS)  fit¬ 
ting  routine  applied  to  the  measured  data.  It  can  be  con¬ 
cluded  from  Fig.  2  that  there  is  a  good  agreement  between 
the  ladder  model  and  the  experimental  data  in  a  wide  range 
of  frequencies.  Also  it  can  be  seen  that  the  imaginary  part 
of  impedance  changes  remarkably  as  the  frequency  varies 
while  the  real  part  keeps  relatively  constant. 

2.7.  Numerical  oscillation 

Typical  solution  techniques  for  differential  equations  in 
most  of  the  time-domain  circuit  simulators  are  trapezoidal 
or  implicit  Euler  or  Gear’s  algorithm,  and  others.  The 
Virtual  Test  Bed  offers  flexibility  for  model  developers  to 
apply  any  of  those  within  their  models  based  on  the  con¬ 
siderations  for  accurate  and  fast  computation.  Since  the 
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Fig.  1.  Equivalent  circuit  of  ultracapacitors. 
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trapezoidal  method  is  absolute  order  stable,  and  provides  a 
good  compromise  among  simplicity  and  accuracy,  it  is  one 
of  the  favorite  methods  used  in  VTB  model  developing. 

However,  numerical  oscillations  can  occur  under  some 
circumstances.  The  numerical  oscillations  can  be  classified 
into  two  types.  The  first  type  of  numerical  oscillation  is 
caused  by  step  change  in  certain  state  variables.  Several  ap¬ 
proaches  have  been  proposed  to  solve  this  type  of  numerical 
oscillation  [9-11].  The  second  type  of  numerical  oscillation 
a-rises  from  an  inappropriate  time  step  size.  Generally,  the 
second  type  of  numerical  oscillation  can  be  avoided  by  re¬ 
ducing  the  simulation  time  step;  however,  unnecessary  small 
simulation  time  steps  will  dramatically  increase  the  compu¬ 
tational  effort.  Clearly,  it  is  necessary  to  choose  an  appro¬ 
priate  time  step  for  better  accuracy  and  faster  computation 
while  avoiding  erroneous  oscillation.  This  paper  focuses  on 
the  second  type  of  numerical  oscillation  in  modeling  and 
simulation  of  the  ultracapacitor  having  multiple  constants, 
and  provides  an  effective  method  to  choose  appropriate  step 
size  for  sufficient  computation  accuracy  while  minimizing 
computation  time. 

Let  us  consider  a  simple  case  first,  a  circuit  including  a 
capacitor  and  a  resistor,  as  shown  in  Fig.  3,  which  may  be 
one  of  the  stages  of  the  ultracapacitor  model. 

The  equation  for  the  terminal  current  is, 


i(t)  =  C 


d  v(t) 
d  t 


d  i(t) 

—  x - 

d  t 


where 


x  =  RC 


is  the  time  constant  of  the  circuit.  Applying  the  trapezoidal 
rule  within  the  time  interval  (t,  t  —  h),  where  h  is  the  time 
step,  and  rearranging,  we  obtain, 

i(t)  =  g-v(t)-b(t-h).  (3) 


where, 


2  C 


g  = 


2r  +  h  ’ 


bit  -h)  = 


2  C  2 r-h 

v(t  —  h) - i(t  —  h) 


2t  T  h 


2x  -f  h 
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(5) 


are  the  admittance  and  the  current  history  respectively.  The 
current  history  is  dependent  of  the  values  of  the  terminal 
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Fig.  3.  A  first  order  equivalent  circuit  of  an  ultracapacitor. 


voltage  and  the  current  at  the  previous  step.  The  following 
recursion  equation  can  be  obtained  using  Eqs.  (3)  and  (5), 


i(nh )  = 


2  x  —  h 
2x  T  h 


m 


v((k  - 1 m  , 


where  nh  =  t.  The  criteria  for  Eq.  (6)  to  be  numerically 
stable  is 

h  >  0,  x  >  0  (7) 


which  can  always  be  satisfied.  However,  if  h  =  2r,  the  first 
term  in  Eq.  (6)  yields  zero.  If  h  >  2r,  the  sign  of  the  first 
term  in  Eq.  (6)  changes  from  one  time  step  to  the  next, 
therefore,  causing  the  current  value  to  oscillate.  The  second 
term  contributes  oscillations  too,  since  the  sign  of  the  voltage 
at  the  Mi  step  also  changes  every  time  n  advances. 

For  the  ultracapacitor  model  shown  in  Fig.  1,  there 
are  total  six  equivalent  time  constants,  corresponding 
to  six  energy- storage  elements,  widely  spreading  from 
micro-seconds  to  seconds.  If  the  model  is  used  in  simula¬ 
tion  of  a  system  of  line-frequency  or  lower  frequencies,  the 
time  step  must  be  smaller  than  two  times  of  smallest  time 
constant  of  the  ultracapacitor,  which  is  in  micro-seconds 
in  order  to  prevent  numerical  oscillations.  For  example, 
Fig.  4  shows  the  numerical  oscillations  of  the  ultracapacitor 
model  with  an  inappropriate  simulation  time  step  (here  us¬ 
ing  0.005  s)  when  it  is  connected  to  an  ideal  voltage  source 
(IV).  It  is  can  be  clearly  seen  that  the  currents  of  C5  and 
C4  oscillate  heavily  from  about  —100  to  100  A  and  decline 
slowly. 

Also  apparent  is  that  the  currents  C\  and  C2  look  much 
better,  which  means  for  a  fixed  circuit  model,  a  given  sim¬ 
ulation  time  step  will  mainly  affect  different  stages.  So  the 
oscillation  problem  can  be  avoided  if  an  equivalent  circuit 
of  an  appropriate  order  or  appropriate  time  constants  is  cho¬ 
sen  based  on  the  characterization  time  of  the  system  inputs, 
as  illustrated  in  the  next  subsections. 


2.2.  Order  reduction  method 

As  shown  in  Fig.  5,  four  different  ladder  networks  are 
embedded  into  one  ultracapacitor  model.  The  model  will 
choose  the  appropriate  circuit  according  to  the  simulation 
time  step.  For  example,  to  reduce  from  the  five-order  to  the 
four-order,  the  4th  order  capacitor  value  in  the  four-order 
network  is  substituted  by  the  C5  +  C4  value  in  the  five-order 
circuit.  The  same  procedure  is  repeated  for  the  resistor  in  the 
four-order  circuit.  Note  the  inductor  is  just  ignored  in  the  low 
frequency  band.  A  close  look  at  the  value  of  capacitor  and 
the  resister  in  the  two-order  shows  they  are  the  summation  of 
all  parameter  values  in  the  five-order.  The  five-order  circuit 
is  not  preferred,  since  it  is  too  simple  and  yields  too  big  error. 
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Fig.  4.  Numerical  oscillation  due  to  an  inappropriate  simulation  time  step. 


Note  that  the  circuit  topology  for  each  order  is  not 
unique  as  it  is  always  somewhat  arbitrary  when  modeling  a 
two- terminal  device  with  more  than  two  internal  elements. 
The  topologies  we  choose,  as  shown  above,  are  simply 
inherited  from  the  original  one,  and  they  are  degraded  to 
a  smaller  RC  network  each  time  an  order  is  reduced.  We 
believe  that  the  network  will  best  represent  the  nature  of 
the  charge  processes  in  the  ultracapacitor,  as  illustrated  by 
the  Fig.  2  and  subsequent  figures. 
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Fig.  5.  Automatic  order  selection  method. 


2.3.  Determination  of  maximum  time  step  for  each  order 


In  this  section,  the  equivalent  circuit  model  in  Fig.  1  is 
first  analyzed.  Then  the  definition  of  the  simulation  time  step 
Tkef  is  proposed;  after  that,  the  model  with  different  order 
networks  as  shown  in  Fig.  5  is  studied  with  consideration  to 
simulation  time  step  selection. 

First,  let  us  consider  the  circuit  shown  in  Fig.  1,  it  is  not 
obvious  how  to  calculate  the  equivalent  time  constants  as¬ 
sociated  with  each  independent  energy  storage  component 
since  all  the  components  are  coupled  together.  However,  the 
time  constants  can  be  approximated  by  decoupling  the  re¬ 
lated  resistor-capacitor  stage  from  the  network  using  an  ap¬ 
propriate  assumption.  Let’s  arbitrarily  choose  node  4  and 
assume  it  is  the  node  whose  variables  first  begin  to  numeri¬ 
cally  oscillate  with  the  simulation  time  step  increasing.  Un¬ 
der  this  condition,  the  voltage  variations  of  the  other  nodes 
are  slower  than  the  voltage  variation  of  node  4.  Approximate 
the  voltages  of  nodes  3  and  5  (V3  and  V5)  as  fixed  values, 
which  means  the  nodes  3  and  5  can  be  looked  upon,  as  con¬ 
necting  to  ideal  voltage  sources,  so  that  the  components  of 
C3,  R2  and  R3  are  decoupled  from  the  network.  The  time 
constant  of  this  stage  resistor-capacitor  is  described  as: 


t3  =  c3 


r3  •  R3 

r3  +  R2 


According  to  Eqs.  (6)  and  (8),  the  maximum  simulation  time 
step  for  node  4  without  numerical  oscillation  is  2T3.  Repeat¬ 
ing  this  process,  all  the  time-constants  of  the  ultracapacitor 
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Table  2 


Six  relevant  time  constants  of  the  model  in  Fig.  1 


Related  components 

Time  constant  (s) 

T6  (R5L) 

9.8187E— 6 

t5  (r5r4c5) 

3.80093E— 5 

T\  (R4R3C4) 

0.001616073 

T3  (R3R2C3) 

0.021411548 

T2  (R2R1C2) 

0.225220521 

Ti  (R1C1) 

0.93832 

model  in  Fig.  1  are  obtained  and  listed  in  Table  2.  From  T& 
to  T\  in  turn,  the  time  constants  increase  respectively,  and 
Tfr  is  the  smallest  one,  which  means  that  the  numerical  os¬ 
cillations  will  occur  from  the  left  stage  to  right  in  turn  as  the 
simulation  time  step  increasing.  This  is  not  a  coincidence 
but  determined  by  the  modeling  method  in  which  the  left 
stages  in  the  model  are  used  to  cover  the  fast  behavior  of  ul¬ 
tracapacitors  and  the  right  stages  are  used  to  determine  the 
long  time  behavior.  21 6  is  the  maximum  time  step  for  the 
model  in  Fig.  1  from  numerical  oscillation;  however,  in  or¬ 
der  to  make  the  simulation  results  represent  the  real  world 
process  well,  the  simulation  time  step  should  be  at  least  a 
few  times  shorter  than  21 6,  we  here  define  the  recommend 
time  step  as: 


2ref  —  5  ^smallest 


where  rsmanest  is  the  smallest  time  constant  in  the  network 
and  Tvq f  is  the  recommend  time  step  for  the  network  simu¬ 
lation. 

After  the  analysis  of  the  equivalent  circuit  model  in  Fig.  1, 
we  consider  the  time  step  selection  for  each  different  order 
network  as  shown  in  Fig.  5.  From  Eq.  (8),  each  time  con¬ 
stant  that  associates  with  each  independent  energy  storage 
component  (capacitor  or  inductor)  in  Fig.  5  can  be  obtained. 
From  these  calculation  results,  we  choose  the  smallest  time 
constants  of  each  different  order  network  and  list  them  in 
Table  3.  Table  3  also  lists  the  recommend  maximum  simu¬ 
lation  time  step  for  each  order  according  to  Eq.  (9). 

It  can  be  seen  from  Table  3,  if  the  five-order  RC  network 
circuit  is  adopted,  the  simulation  time  step  should  be  set 
smaller  than  7j-ef5 ;  if  the  four-order  network  circuit  is  used, 
the  simulation  time  step  should  be  limited  to  the  range  be¬ 
tween  rref5  and  rref4,  etc.  The  maximum  simulation  time 
step  for  this  model  is  7j-ef2.  Notice  that  if  a  two-order  is  cho¬ 
sen  for  simulation,  the  time  step  size  now  is  about  0.16  s, 


Frequency  (Hz) 

Fig.  6.  Comparison  of  different  order  model  simulation  results  ( — )  with 
measurement  (□  and  *). 

while  the  original  model  needs  a  step  size  of  2  p,s.  Signifi¬ 
cant  computational  time  saving  is  obvious. 

2.4.  Model  order  selection 

Fig.  6  shows  the  comparison  among  different  order 
networks  simulation  results  and  the  measurements.  The 
five-order  fits  the  whole  test  frequency  spectrum.  The  imag¬ 
inary  part  the  complex  resistance  first  decreases  as  the 
frequency  increases,  and  then  it  increases  because  of  the 
series  inductor.  The  four-,  three-  and  two-order  fit  the  low 
frequency  bands  well  respectively. 

To  discuss  the  merits  of  the  reduced  order  model,  we  cal¬ 
culate  the  complex  impedance  error  between  two  adjacent 
orders.  As  shown  in  Fig.  6,  the  error  between  two  consecu¬ 
tive  orders  increases  as  the  frequency  increases.  Define  this 
error  as: 

771  /  s\  I I  r^jri  1  ( f)  I  I  '^jn  ( /)  1 1  |  f\r\o7  /  t 

Enxn-l){f)  =  - .  ,A| -  X  100%  (10) 

I  \J)  I 

where  En^n-\)(j)  is  the  error  between  order  n  and  order 
n  —  1  at  frequency/,  n  is  the  order  number,  Zn  is  the  complex 
impedance  and /is  the  frequency. 

The  frequency  at  which  the  error  exceeds  1  %  can  be  cal¬ 
culated  according  to  Eq.  (10)  and  they  are  listed  in  Table  4. 

Table  5  combines  the  Tables  3  and  4,  which  explains  the 
model  order  selection  with  maximum  1%  truncation  error. 


Table  3 


The  smallest  time  constant  in 

each  order  network  and  the  recommend  time  steps 

Different  order  network 

Related  components  of  the  smallest  time  constant 

Time  constant  (s) 

Recommend  maximum 
simulation  time  step  (s) 

Five-order  RC  +  L 

Four- order  RC 

Three-order  RC 

Two-order  RC 

R5  -  L 

(R5  +  R4)  —  (C5  +  C4)  —  R3 

(P5  +  R4  AR3 )  —  (C5  +  C4  +  C3)  —R2 

(R5  R4  R3  +  Rf)  ~~  (C5  +  C4  +  C3  +  Cf)  —R\ 

9.8187E— 6 

0.002596 

0.0688768 

0.729579 

rref5  =  1.964E— 6 
rref4  =  5.192E— 4 
rref3  =  0.01378 
rref2  =  0.145916 
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Table  4 


Frequency  at  which  error  >1% 

Error 

1%  Error  frequency  (Hz) 

E5-4 

550 

£4-3 

1.44 

E3-2 

0.0034 

Table  5 

Model  order  selection  method 

Frequency  band  of 

Model  order 

Recommended 

interest  (Hz) 

simulation  time  step  (s) 

/>  550 

6 

T  <  1.964E— 6 

1.44  </<  550 

4 

1.964E— 6  <  T  <  5.192E— 4 

0.0034  </<  1.44 

3 

5.192E— 4  <  T  <  0.01378 

/  <  0.0034 

2 

0.01378  <  T  <  0.145916 

For  example,  if  the  relevant  frequency  band  of  a  simulation 
is  wider  than  550  Hz,  the  ultracapacitor  model  will  adopt  the 
five-order  RC  circuit,  and  the  recommend  time  step  for  this 
simulation  is  T  <  1.964E— 6  s. 

Fig.  7  shows  the  working  frequency  band  for  the  different 
order  networks. 


2.5.  Capacity  scaling 


Another  useful  characteristic  of  this  model  is  capacity 
scaling.  Model  default  parameter  values  are  based  on  a 
Maxwell  100  F  ultracapacitor,  but  the  model  will  automati¬ 
cally  calculate  new  parameters  according  to  users’  specified 
capacitance  by  using  Eq.  (11): 


N  = 


CTotal 


100  ’ 

Ck_ NEW  =  Ck  •  N, 


Rk- NEW  = 


Rk 

N 


L 

^NEW  =  — 


(ID 


where  N  is  the  capacity- scaling  ratio,  and  74_new,  Q  _new 
and  Lnew  are  the  new  values  corresponding  to  the  new 
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Fig.  7.  Working  frequency  band  of  different  orders. 


total  capacitance.  This  scaling  is  a  little  risky  considering 
that  electrode  parameters  may  change  with  capacity,  but  it 
nonetheless  gives  a  quite  good  estimate  of  capacitor  perfor¬ 
mance  for  system-level  studies.  Fig.  8  shows  the  accuracy 
of  scaling  the  capacity  from  100  F  (default  value)  to  10  and 
5  F,  and  comparison  with  measurements  of  Maxwell  10  and 
5  F  capacitors.  On  the  left  is  the  comparison  of  Maxwell 
10  F  and  on  the  right  is  5  F. 

From  Fig.  8,  we  can  see  that  the  model  with  capacity 
scaling  matches  well  experiments  within  the  same  brand 
product  family.  Comparison  with  different  brand  products 
may  have  a  larger  error. 


3.  Model  used  in  FCHV  simulation 

One  important  application  of  ultracapacitors  is  in  the  fuel 
cell  powered  hybrid  electric  vehicles  (FCHVs),  in  which 
the  ultracapacitor  provides  peak  power  during  acceleration 
or  serves  as  the  immediate  repository  for  electric  energy 
regenerated  during  braking.  Simulation  of  a  hybrid  electric 
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DC  Bus 


*D.C.—  Distributed  Controller 


Fig.  9.  Block  diagram  of  hybrid  electric  vehicle. 

vehicle  is  a  complex  process  involving  many  technical  disci¬ 
plines.  The  interdisciplinary  nature  of  the  problem  requires 
its  component  models  to  be  accurate,  fast  and  robust. 

The  FCHV  described  in  this  paper  is  composed  of  a  me¬ 
chanical  train  that  includes  all  the  mechanical  parts  of  the 
vehicle,  including  the  motor  that  converts  energy  between 
electrical  and  mechanical,  a  fuel  cell  system  as  the  primary 
source  of  power,  battery  and  ultracapacitor  stacks  to  meet 
high  and  intense  power  demands  respectively,  D.C./D.C. 
converters  to  control  power  flow  between  the  components, 
and  a  supervisory  controller  to  control  and  balance  the 
whole  system.  Fig.  9  shows  the  block  diagram  of  the 
FCHV. 


Fig.  11.  Speed  of  vehicle  during  EPA  city  driving  test. 

As  shown  in  Fig.  9,  the  ultracapacitor  bank  is  connected 
to  the  D.C.  bus  through  a  bi-directional  D.C./D.C.  converter, 
which  is  controlled  by  the  vehicle  supervisory  controller. 
Fig.  10  shows  the  FCHV  system  schematic  in  the  VTB  en¬ 
vironment. 

The  designed  FCHV  is  tested  using  the  Environmental 
Protection  Agency  (EPA)  city  driving  test  profile.  Some  sim¬ 
ulation  results  are  shown  in  Figs.  11  and  12.  Fig.  11  shows 


Fuel  Cell  Sys. 


Chassis  Dynamics 


fOMDl 

poo  — «M> 
pooj 

EPA  Cycle 


Fig.  10.  Schematic  of  FCHV  system  in  VTB. 
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Fig.  12.  Terminal  voltage  of  the  ultracapacitor  bank  during  EPA  city 
driving  test. 

the  vehicle  speed  variation  during  the  EPA  city  driving  test 
with  the  y-axis  as  the  speed  in  mil/h  and  the  v-axis  as  the 
time  in  second.  It  can  be  seen  the  vehicle  speed  (dash  curve) 
follows  the  commands  (solid  curve)  very  well. 

Fig.  12  shows  the  simulation  result  of  the  ultracapacitor 
voltage  variation  during  the  EPA  city  driving  test,  the  y-axis 
is  the  voltage  in  volts  and  v-axis  is  the  time  in  seconds.  The 
total  capacity  of  the  ultracapacitor  model  is  scaled  from  100 
to  76  F  according  to  the  vehicle  energy  requirement.  The 
initial  value  of  the  ultracapacitor  terminal  voltage  is  set  as 
600  V  in  order  to  store  more  energy  and  in  turn  to  decrease 
the  usage  of  ultracapacitors.  The  ultracapacitor  voltage  in¬ 
creases  to  the  maximum  value  about  620  V  when  the  vehi¬ 
cle  decelerates,  and  decreases  to  about  the  minimum  value 
about  540  V  when  the  vehicle  accelerates. 

In  this  simulation,  the  simulation  time  step  is  set  as  0.01  s. 
According  to  Table  5,  the  ultracapacitor  works  in  three-order. 
The  simulation  works  well  without  numerical  oscillations. 


4.  Conclusions 

A  multi-stage  resistantce/capacitor  ladder  model  has  been 
built  with  automatic  order  selection  and  capacity  scaling 
for  ultracapacitor  modeling.  The  feature  of  the  automatic 
order  selection  makes  the  model  automatically  select  the 
order  or  number  of  stages  of  the  ladder  network  according 
to  the  simulation  time  step  to  minimize  calculation  error 


and  save  computational  cost.  The  feature  of  the  capacity 
scaling  makes  the  model  accurately  reflect  other  capacitors 
with  different  sizes.  A  model  with  automatic  order  changing 
up  and  down  while  the  simulation  is  running  with  variable 
simulation  time  step,  will  be  explored  in  future  work. 
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